###########################CARREGAR PACOTES##############################
#########################################################################
require(RcmdrMisc) #pacote para carregar arquivos do Excel
require(agricolae)
require(dae)
require(ExpDes)


#Mudar diretório
setwd("D:/OneDrive/UFMT/Disciplinas/Veterinaria/experimental/2019/dados")
setwd("C:/Users/ander/OneDrive/UFMT/Disciplinas/Veterinaria/experimental/2019/dados")
source("funcoes.txt")



#########################################################################
############################Análise do DBC###############################
#########################################################################

##Importar dados
dados=readXL("coelhos.xls",sheet = 1,stringsAsFactors=TRUE)
attach(dados)
str(dados)


Bloco=as.factor(rep)
Trat=as.factor(sub)
VR=peso
######Análise de Variância#
model=aov(VR~Bloco+Trat)
anova(model)
cv.model(model)

#####PRESUPOSIÇÕES DA ANOVA
#Gerar os erros padronizados do modelo
e=rstandard(model)



#Teste de Shapiro-Wilk - NORMALIDADE
shapiro.test(e) 
#Teste de Oneill-Mathews - Homogeneidade de Variância
oneillmathewsDBC(VR,Trat,Bloco)

##Teste de Durbin-Watson - Independencia
set.seed(12345)
durbinWatsonTest(model)

###Teste de aditividade dos efeitos de tratamento e bloco
df<-df.residual(model)
MSerror<-deviance(model)/df
nonadditivity(VR, Bloco, Trat, df, MSerror)


####Anova com Teste de Tukey####################
rbd(sub,rep, peso, quali = FALSE, sigT = 0.05, sigF = 0.05)

grafico.reg(sub,peso,grau=1,xlab="Substituição de Proteina",ylab="Peso")